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Abstract. We develop a method to measure the strength of the absorption features in Type la supernova (SN la) 
spectra and use it to make a quantitative comparison between the spectra of Type la supernovae at low and high 
redshifts. In this case study, we apply the method to 12 high-redshift (0.212 < z < 0.912) SNe la observed by the 
Supernova Cosmology Project . Through measurements of the strengths of these features and of the blueshift of 
the absorption minimum in Ca II H&K, we show that the spectra of the high-redshift SNe la are quantitatively 
similar to spectra of nearby SNe la (z < 0.15). One supernova in our high redshift sample, SN 2002fd at z=0.279, 
is found to have spectral characteristics that are associated with peculiar SN 1991T/SN 1999aa-like supernovae. 



1. Introduction 

Type la supernovae (SNe la) are excellent distance indica- 
tors and have been used to show that th e expansion of the 
Universe is currently acceler at ing (iPerlmutter et al 



19981 iGarnavich et all 1 19981 I Schmidt et alj 1 19981 




Perlmutter et al.l 19991: iKnop et al 



Krisciunas et al.l 
2007). 
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Wood-Vasev et al 



2005; 



Astier et al 



2006; 
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Recently, both the SNLS0 and ESSENCE^ projects, 
which aim to use hundreds of SNe la to constrain the 
nature of dark energy through the measurement of the 
equation-of-state parameter, have reported their first re- 
sults. The control of the systematic uncertainties is criti- 
cal to their success. Among the possible systematic effects, 
evolution of the SNe la population over cosmological time- 
scales is one of the most important and least understood. 

Spectra of SNe la are well-suited to study potential 
evolutionary effects. For example, the average metallicity 
of the Universe increases with cosmic time, so it is reason- 
able to expect that high-redshift SNe la are in environ- 
ments that have lower average metallicity than those of 
nearby SNe la. The effect on the spectral energy distribu- 
tion of a lower metallicity progenitor has bee n modeled by 
Hoefiich et ail (|l998l ) and lLentz et~ai1 (|2000h . These stud- 
ies find that such SNe la, especially at early epochs, are 
expected to show enhanced flux in the UV, weaker absorp- 
tion features in the optical and a shift in the absorption 
minima of optical features to longer wavelengths. 

With the large number of well-observed low-redshift 
supernovae now availab le, a wide ra n ge of spectral di- 



versity is being found (jBranch et al. ( 20061 see e.g. )) 



The physical origin of these differences is still not com- 
pletely understood making it difficult to predict their pos- 
sible evolution with redshift. Statistical studies are use- 
ful to probe differences between high and low-redshift 
SN la data sets. So far, few distant SN la spectra have 
been compared with low-redshift data sets in a quanti- 



Howell et al.l 



1998; Coil et al 



200C 



Barris et al. 2004 Riess et al. 2003 


Blakeslee et al.ll2003t 


Matheson et al.1 2005: Hook et al 


12005: iBalland et al 



Blondin ct al 



2006), 



and few 

spectroscopically confirmed high-redshift SNe la have 
been repo rted as peculiar ( two S N 1991T/SN 1999aa-like 



SNe la in Matheson et al. (2005); one, as we will see 



the present paper). Hereafter, we follow the convention 
of describing SN 1991T/SN 1999aa-like SNe la by using 
"91T-like" to represent both SN la subtypes. 

During 2000, 2001 and 2002, the Supernova Cosmology 
Project (SCP) obtained spectra of 20 high-redshift 
SNe la with FORS2 on the ESO Very Large Telescope 
( Lidman et all 120051 ). [n this paper, we analyze the 14 



spectra of 12 SN la with the highest signal-to- noise ra- 
tios and perform a quantitative comparison between these 
spectra and the spectra of low-redshift SNe la. The def- 
inition of a newly introduced spectral indicator for Type 
la supernovae, namely pseudo-equivalent width, is given 
in section [5] The data-sets we apply the method to are 
presented in section [3l The properties of the pseudo- 
equivalent width for SNe la are given in section[4]together 
with the result of the comparison between high- and low- 
redshift SNe la. Summary and conclusions are given in 
section [5] 



2. Method 

We use two spectral indicators to compare high- and low- 
redshift SN la spectra: 

1) the wavelength of the absorption minimum at ~ 
3750 A, which is largely due to Ca n H&K, and 

2) newly-introduced pseudo-EW measurements of fea- 
tures associated with Ca n H&K, Mg II and Fe II. 

2.1. Pseudo-Equivalent Widths of SNe la 

The spectra of Type la supernovae are very characteris- 
tic and, in comparison with supernovae of other types, 
relatively homogeneous. Spectral features are broad, re- 
flecting the high velocities of the ejecta (~10000 km s _1 ), 
and evolve with phase. However, differences b etween dif- 
ferent Type la supernovae have been noted (see Filippenkol 
(jl997l ) for a review) . In some cases the differences are dra- 
matic and have resulted in the definition of SNe la sub- 
types. 

The p seudo-equivalen t width (pseudo-Ew), first de- 
scribed in iFolatellil (|2004l ). can be used as a spectral indi- 
cator. It is described in detail in the following section. 

2.2. Definition 

The equivalent width, as used in stellar spectroscopy, can 
be used to measure the strength of absorption features in 
supernova spectra. However, the relationship between this 
quantity and the physical properties of the SN la ejecta 
is complex. In a SN la spectrum the overlap of thousands 
of lines give rise to a "pseudo" continuum thus a real con- 
tinuum can not be identified. Because of this distinction 
with the definition of equivalent widths we use the term 
"pseudo-equivalent widths " to refer to our measurements. 
This does not prevent us from using this well-defined quan- 
tity, nor does it not allow up from building a consistent 
set of measurements for all SNe la in our data-sets. 

In this analysis, we define an absorption feature as a 
wavelength region that is bounded by two local flux max- 
ima. Fig.[Tjshows eight such regions, corresponding to the 
eight strongest absorption features at optical wavelengths 
in supernova spectra up to approximately one month af- 
ter maximum light. Each feature is marked with a number 
from 1 to 8 and each number corresponds to a mnemonic 
name. The upper and lower limits of the features vary 
in time (because of the SN la envelope expansion), and 
from SN la to SN la at a given phase (i.e. SN la spectral 
diversity). We give ranges for these limits in Table [T] 

As already mentioned, in order to measure a pseudo- 
ew, a pseudo-continuum must be determined. We define 
the pseudo-continuum as the straight line fit through the 
two local maxima that bound a feature^ (see panel (b) in 
Fig. []}. A detailed explanation of the measurement tech- 
nique is given in section l2~3l 



1 http://www.cfht.hawaii.edu/SNLS/ 

2 http:/ /www. ctio.noao.edu/essence/ 



The chosen maxima are those that maximize the wave- 
length span of the feature with the restriction that the derived 
pseudo-continuum does not intersect the spectrum within the 
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Once the pseudo-continuum is defined, the pseudo-EW 
is computed for each feature within its wavelength limits. 
The spectrum is divided by the pseudo-continuum and the 
resulting area of the feature is measured (in units of A). 
In this case, the calculation was approximated by a simple 
rectangular integration method: 



pseudo-EW = 2_. ( 1 



fx (AQ 
/c(Ai) 



AA,, 



(i) ¥ 



where Xi (i = 1, . . . , N) are the central wavelengths of the 
bins of size AAi, f\(Xi) is the measured flux in bin i, and 
fc(Xi) is the fitted pseudo-continuum flux evaluated at the 
same points. Deviating points due to bad pixels, night sky 
residuals, or narrow host-galaxy lines were rejected using 
a 3<7-clipping algorithm. 




MA) 



8000 9000 



Table 1 Feature limits. 



Feature 
ID 



Mnemonic 
Label 



Blue- ward 
limit range (A) 



Red-ward 
limit range (A) 



1 


"Ca ii H&K" 


3500 


- 3800 


3900 


- 4100 


2 


"Si ii 4000" 


3900 


-4000 


4000 


- 4150 


3 


"Mg ii 4300" 


3900 


- 4150 


4450 


- 4700 


4 


"Fe ii 4800" 


4500 


- 4700 


5050 


- 5550 


5 


"S ii W" 


5150 


- 5300 


5500 


- 5700 


6 


"Si ii 5800" 


5550 


- 5700 


5800 


- 6000 


7 


"Si ii 6150" 


5800 


- 6000 


6200 


- 6600 


8 


"Ca ii IR" 


7500 


- 8000 


8200 


- 8900 



The la statistical uncertainty was computed by error 
propagation from the estimated uncertainties in the spec- 
tral flux (erf) and in the pseudo-continuum (cr c ): 

■ 



^pscudo-EW 



■ N 



/a 2 (A.) , 



8=1 



/ C 2 (A*) / C 4 (A 4 ) 



(A,) (AXiY 



(2) 

By definition, the rest-frame pseudo-equivalent width 
is redshift independent and thus provides a useful tool for 
comparing spectra from different objects over a wide range 
of redshifts. 



2.3. Measurement technique 

Measuring pseudo-equivalent widths on high signal-to- 
noise ratio spectra is relatively simple since the local 
maxima bounding an absorption feature can be identi- 
fied easily. On low signal-to-noise data the measurements 
are more difficult. For comparing high redshift supernovae 
(which generally have low signal-to-noise) with local SN la 
spectra, we have established a measurement technique (to 
be applied on spectra regardless of their signal-to-noise 
ratio) that minimizes possible systematic effects. 

To measure the pseudo-EW of a spectral feature the 
local pseudo-continuum must be determined. To perform 
this operation we proceeded as follows: 



feature limits with the possible exception of noise artifacts in 
case of low signal-to-noise ratio data. 



Fig. 1 SNe la spectral evolution and feature definitions 
for three epochs: 2, 16 and 39 days after maximum light. 
Numerical labels correspond to the following adopted fea- 
ture names: 1- "Ca II H&K"; 2- "Si n 4000"; 3- "Mg II 
4300"; 4- "Fe n 4800"; 5- "S II W"; 6- "Si II 5800"; 7- "Si II 
6150"; and 8- "Ca II IR". Short vertical lines show the ap- 
proximate positions where the pseudo-continuum is taken 
in each case. Feature ranges change with time and, due to 
blending, some weaker features are not considered at later 
epochs. Note that after 2-3 weeks past maximum light, the 
selected pseudo-continuum points correspond to emission 
peaks. Panel (a): the region around features #2 and #3 
for near-maximum spectra of SN 1991T (top), SN 1989B 
(middle), and SN 1991bg (bottom). Feature #2 is not de- 
fined in the case of 1991bg-like SNe la because the region is 
dominated by absorption from Till. Adopted feature lim- 
its are marked with vertical lines. Panel (b): an example 
of the pseudo-continuum trace for "Fe II 4800" on a normal 
SN la near the time of maximum light. Here, solid verti- 
cal lines show the regions where the pseudo-continuum is 
fitted. Dotted lines mark the bounds used to measure the 
pseudo-EW. 



— The two local maxima that bound the absorption fea- 
ture are visually identified and marked. These are the 
bounds used in performing the sum in Eqn. [1] 

— A small wavelength region (to which hereafter we refer 
to as the fitting region) is selected around each identi- 
fied maximuny, always within the wavelength ranges 
listed in Table [Q 

— A straight line is fitted to the data. The result of the 
fit is taken to be the pseudo-continuum and is used in 
Eqn. [Q 

The wavelength span of the two fitting regions depends 
on the morphology of the local maxima. On the blue end of 
the spectrum, where more absorption features are found, 
the maxima tend to be narrower than on the red end. 



4 See panel (b) in Fig. [T] the four vertical lines represent an 
example of the wavelength span of each fitting region. 



4 



G. Garavini: Search for Spectral Evolution in high-redshift supernovae 



Thus, the fitting regions are, on average, smaller (~ 20 A) 
in the blue and slightly larger (~ 25 A) in the red part 
of the spectrum. Larger wavelength spans would tend to 
set a lower pseudo-continuum level and thus to produce 
pseudo-EW systematically biased toward smaller values. 
The size of the fitting region is independent of the signal- 
to-noise ratio; however, in low signal-to-noise ratio spec- 
tra, noise spikes that could bias the pseudo-continuum fit 
are clipped. 

3. Data sets 

In the following section, we present the high and low red- 
shift data-sets that we have used in this study. 



3.1. High Red shift Data-set 

The supernova spectra that are analyzed in this work were 
obtained as part of several campaigns by the SCP to dis- 
cover and follow a lar ge number of S Ne la over a wide 
range of redshifts (see ILidman et al] ( 2005) for details) 



Out o f the 20 spectrally confirmed SNe la in ILidman et al 
( 20051 ). we select the 12 SNe la (z=0. 212-0.912) with 



signal-to-noise ratios S/N>3 per 20 A bin to pursue our 
quantitative analysis. One supernova, SN 2001go, was ob- 
served at three epochs, so there are 14 spectra in total. 

The biases affecting this sample are complex. First, 
supernova searches are magnitude limited, so Malmquist 
and light curve shape biases (i.e. brighter/dimmer SN la 
have broader/narrower light curves) make it unlikely that 
low-luminosity, SN 1991bg-like SNe la will be found. Also, 
it is likely that, while ranking supernovae candidates for 
follow-up spectroscopy, 1991bg-like SNe la would be as- 
signed lower priority due to their low luminosity. Second, it 
is more difficult to spectrally confirm high-redshift SNe la 
in spectra where host galaxy contamin ation is above 75% 
( Lidman et alll2005l : iHowell et al.ll2005l ). Although this af- 
fects all SNe la to some degree, lower-luminosity SNe la 
are more commonly found in bright ellipticals Moreover, 
we can not exclude that the signal-to-noise cut we per- 
formed while choosing the sub-sample of spectra to study 
did not introduce an extra selection bias toward over- 
luminous objects. As it will be clear in sectional the spec- 
tral indicators we use to search for evolution with red- 
shift do not unambiguously discern between normal and 
over-luminous SN la. The additional scatter that would 
be introduced by using data with lower signal-to-noise ra- 
tios would obfuscate the result. We make no attempt to 
correct the sample for these biases. 

The spectra, re-binned to 20 A are shown in Fig. [21 and 
are summarized in Table [21 A full descriptio n of the obser 



yation s and the data reduction are given in ILidman et al 



(2003). For each supernova, we use the error spectrum, 
which was estimated from regions free of SN la and host 
galaxy light on the sky-subtracted two dimensional spec- 
trum, to estimate the statistical uncertainties on the quan- 
tities we compute in the following analysis. 



High-redshift supernova spectra can contain significant 
amounts of light from the host galaxy. On the 2d spectra, 
the host galaxy and the SN la are often spatially unre- 
solved, making it difficult to estimate the contribution of 
the host to the observed flux. We estimated this contri- 
bution using a template matching technique based on a 
large set of nearby supe rnovae spectra and g alaxy models 



similar to those used in Lid man et al.l ( 20051 ). The contri- 



bution of the galaxy light, relative to the total observed 
flux, are tabulated in column 10 of Table [21 

The epochs with respect to the B— band maximum 
light were estimated using the preliminary light curves, if 
available, and/or spectro scopic dating by template match- 



ing w ith low-z SNe la (jLidman et al.l 120051 : iHook et alJ 
200 51). The two methods usually agr ee within three days 



(jHook et al]|2005l : lHowell et al.ll2005l) . therefore we take 3 
days to be the uncertainty on the quoted epoch whenever 
a light curve estimate of the maximum was not available. 
The redshift of the supernova, when quoted with 3 signif- 
icant figures, was estimated from host galaxy lines visible 
in the spectrum. When this was not possible, the redshift 
was estimated from supernova spectral features, and is 
then quoted with 2 significant figures to account for the 
large intrinsic width of SN la spectral features. 

3.2. SN la identification 

The identification of SN la relies primarily on the detec- 
tion of the absorption feature at approximately 6150 A 
due to Si II A6355. At redshifts above z=0.5, however, this 
characteristic feature is redshifted beyond the wavelength 
range of most optical spectrographs and the classification 
of the supernova ha s to rely on spectral features that lie at 
bluer wavelengths (ILidman et al. 20051 : Hook et al. 2005 : 
Matheson et a l. 2005). Because of the low signal-to-noise 
ratio usually available in high-redshift supernova spectra, 
this approach is not always conclusive. 

We can also use the spectra to identify spectral pe- 
culiarities among SNe la as those found in 91T-like or 
91bg-like supernova. In TableEl the characteristics of four 
wavelength regions for different types and sub- types of su- 
pernovae are schematically reported. Each spectral feature 
is qualitatively described as strong, weak or absent based 
on the absorption strength and broad or narrow based on 
the wavelength span. In the absence of a procedure that 
is based on quantitative measurements, this scheme helps 
in identifying the SN type and, in the case of SNe la, the 
sub- type. 

3.3. SN 2002 fd: A SN 1991T/SN 1999a a-l ike 
Supernova 

SN 2002fd (z = 0.279) is the only supernova in our 
data set that clearly deviates from a "normal la" (see 
Table [3]) . Within the scheme described above, the spec- 
trum of SN 2002fd is similar to the spectra of SN 1999aa 
(jGaravini et al.ll2004l ). a 91T-like supernovae (Fig[3]). The 
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I 1 1 , , 1 1 1 1 

2000 3000 4000 5000 6000 7000 8000 

Wavelength [A] 

Fig. 2 Spectra of the high-redshift SNe la (re-binned to 20 A ) used in this study, plotted in rest frame. For each 
spectrum, we indicate the redshift and epoch (in days from B-band maximum; square brackets). See Table[5]for details. 
Gray vertical bands approximately indicate the wavelength regions used for the quantitative comparison presented in 
section S) 
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Table 2 A summary of the high-redshift data. (See text for details.) 



SN- 


name 


Redshift a 


Date 
(MJD) 


Days from b 

B-band 
Maximum 


Instrument 




setup 




telescope 


Exposure 
time (s) 


S/N c 


Galaxy d 

% 


SN 


2001gy 


0.511 


52021.3 


-7,5(1) 


FORS1 


300V 


grism 


+ 


GG435 


VLT-UT1 


2400 


11 


20 


SN 


2002fd ° 


0.279 


52376.1 


-7(3) 


FORS2 


300V 


grism 


+ 


GG435 


VLT-UT4 


600 


46 


28 


SN 


2000fr 


0.543 


51676.2 


-5.5(1) 


FORS1 


300V 


grism 


+ 


GG435 


TT m T Tm-i 

VbT-UTl 


7200 


21 


16 


SN 


2001gw 


0.363 


52021.4 


-1(3) 


FORS1 


300V 


grism 


+ 


GG435 


TT m T Tmi 

VbT-UTl 


1200 


8 


19 


SN 


2001ha 


0.58 


52022.0 


0(3) 


FORS1 


300V 


grism 


+ 


GG435 


t TT TTmi 

VLT-UTl 


3600 


4 


6 


SN 


2002gl 


0.510 


52413.1 


0(3) 


FORS2 


300V 


grism 


+ 


GG435 


t TT m T Tm A 

VLT-UT4 


3000 


9 


23 




zuuinc 


U.oo 


conoo 1 
OZUZZ. 1 


0(3) 


r Uilol 


oUU V 


grism 


+ 




\n T TTT 1 ! 
V Ll - U 11 


loUU 


1 A 
14 


14 


SN 


2001gr 


0.540 


52021.0 


2(3) 


FORSl 


300V 


grism 


+ 


GG435 


VbT-UTl 


3600 


5 


57 


SN 


2002gi 


0.912 


52407.2 


2(3) 


FORS1 


3001 


grism 


+ 


OG590 


VFT-UT3 


7200 


3 


37 


SN 


2001gm 


0.478 


52021.3 


5(3) 


FORSl 


300V 


grism 


+ 


GG435 


VbT-UTl 


2400 


3 


28 


SN 


2001go 


0.552 


52021.3 


5.6(1) 


FORSl 


300V 


grism 


+ 


GG435 


VbT-UTl 


2400 


5 


13 


SN 


2002gk 


0.212 


52413.3 


6(3) 


FORS2 


300V 


grism 


+ 


GG435 


VbT-UT4 


900 


34 


66 


SN 


2001go 


0.552 


52027.1 


9.5(1) 


FORSl 


300V 


grism 


+ 


GG435 


VbT-UTl 


7200 


10 


19 


SN 


2001go 


0.552 


52058.1 


29.5(1) 


FORSl 


300V 


grism 


+ 


GG435 


VbT-UTl 


9000 


4 


53 



a The redshift is quoted to three significant figures if it is determined from host galaxy lines; two significant figures when determined 
from SN la spectral features. 

''Uncertainties are quoted in parentheses. When the epoch is quoted with 3 days uncertainty it refers to that of the best matching 
low-z la template. 
c per 20 A bin 

d Estimated fraction of the galaxy light in the SN la spectrum, expressed as percentage of the total flux. 
e 9 IT- like Type la SN 

Table 3 A description of the spectroscopic features used to type SN la at z > 0.5. Four wavelength regions are selected 
for performing the SN la typing. See the text for details. 

Region Id A-Region Normal Type la 91T/99aa-like 91bg/86G-like Type Ib/c 

Rest Frame [A] 

'Ca n H&K' a 3550-3950 strong/broad weak or absent/broad b,c strong/broad evident /broad cl 

'Si ll' c 3950-4100 evident/narrow c weak b,c absent absent 

'Fe ll' f ' 9 4600-5200 strong/broad strong/narrow 15 strong/broad strong/broad 

'S II W 5200-5600 strong/narrow b,c weak or absent b strong/narrow b,c absent 



"High and low velocity component in some SNe la. 
6 Around 1 week before maximum light 
c Around maximum light 
d A few exceptions. 

e Marks the beginning of the distinctive strong Ti II absorption feature in 91bg-like SNe la. 
^In Normal la characteristic line profile time evolution. 
s Dominated by Fe ill in pre-maximum spectra of 91T-like SNe la. 



strength of the 'Ca n H&K' feature in SN 2002fd is 
stronger than in 91T-like SNe la and weaker than in nor- 
mal SNe la. The 'Fe ll' and 'S II W reg i ons a re simi- 



lar to those in SN 1999ac (|Garavini et al.ll2005T ). Given 
the low-redshift, Si II A 6355 is also visible. This feature 
in SN 2002fd is intermediate in strength compared with 
SN 1999aa and in normal SNe la. From this qualitative 
analysis of the spectrum, we classify SN 2002fd as a pecu- 
liar SN la, similar to the 91T-like SN 1999aa. In section 



14.31 we show that the pseudo-equivalent widths of the ab- 
sorption features in SN 2002fd are also consistent with 
those found at low redshift for the 91T-like objects. 

Finding SNe la with spectral characteristics similar 
to those of SN 1991T at high redshifts is important. 
The lack of such SNe in high redshift surveys might 
be a sign of e volution. In a distance limited survey, 
iLi et al.l (|2001bh found that approximately 20% of the 
analyzed data set of nearby SNe la could be classified 
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SN 2002fd -7 




4000 



5000 6000 
Rest Frame X[A] 



7000 



Fig. 3 SN 2002fd at day -7, re-binned to 20 A per pixel 
(thick solid line), compared to normal and peculiar SNe la. 
The © symbols mark regions of strong telluric absorption. 



as SN 1991T/SN 1999aa-like which w ere peculiar, over 



luminous SNe la ( Garavini et al. 2004). This percentage 



is probably higher than that of peculiar 91T-like found 
at high-redshift. Cur rently, the reported fraction of 91T 



like S Ne la is < 5% (jLidman et all 1200a iMatheson et al. 
l2005h . However, the fraction of these SNe la in high- 
redshift surveys is uncertain because of the difficulty in 
identifying such SNe la. Several more 91T-like SNe la 
may have already been observed at high-redshift but not 
clearly identified because the spectrum was taken well af- 
ter maximum light or the light-curve had a Amis consis- 
tent with normal SNe la. Over-luminous SNe la such as 
91T-like SNe la generally have broader light-curves than 
normal SNe la. However, it is now becoming evident that 
91T-like SNe la do not always have high Amis values 
and that broad light-curve SNe la do not always show 
the spectra peculiarities seen in 91T- like SNe la. Exampl e 



of the latte r case are SN 1999ee (lHamuy et al.l l2002h 



SN 2 02CX dLi et al.ll2003h or SN 1999aw istrolger et al'l 
20021 ). Li et al.l ( 2001a ) computed that the peculiarity rate 
for 91T-like SNe la should vary between 6% and 18.6%. 
In this case th e observed pecul i arity r ate at high redshift 



(about 4% in IMatheson et all (|2005h and about 5% in 



liar supernovae is highly dependent on the search strategy. 
Nevertheless, the identification of peculiar SNe la in high- 
redshift samples is an important step toward determining 
whether the range of SNe la sub-types that is observed at 
low redshift is also observed at high redshifts. Of course, 
large data samples with high signal-to-noise ratio spectra 
are required to check whether the fraction of peculiar high- 
redshift SNe la is consistent with that found in the low- 
redshift Universe or if there is an evolution with redshift 
in the relative fractions, which might affect the derivation 
of cosmological parameters from SNe la. 

3.4. Nearby Supernova Data-sets 

Two different data sets are used to derive the properties 
of pseudo-EW in local SNe la. 
Set A 



77 spectra (presented in lFolatellil (|2004h ) from 
13 of the SNe la discovered and followed by the Supernova 
Cosmology P roject (scp) in co llaboration with members 



of th e eros ([Hardin et al.l 120001 ) . quest (ISchaefer et al 
1999), and Nearby Galaxies SN Search (|Gal-Yam et al 



our data set) is consistent with the computations, espe- 
cially if one considers that the expected number of pecu- 



1999) teams. The redshift range of these SNe la is 0.01 < 
z < 0.15. 

The two-dimensional raw images were reduced according 
to standard procedures. The calibrated spectra were ad- 
ditio nally corrected for atmospheric and Galac tic extinc- 
tion (jCardelli et al.lll989l : ISchlegel et aD 1 19981 ) and their 
flux-calibration was checked against measured broad-band 
photometry and found in agreement within the quoted un- 
certainties. The spectra were de-redshifted. More details 
about these data can be found in Table [4j All the spectra 
in Set A include estimated statistical uncertainties for 
each wavelength bin. Host-galaxy light was present in a 
subset of these spectra. This contribution was estimated 
and subtracted (for details on the procedure see Sec. 13. ip 
in those cases where it exceeded 10% of the total flux. 

Set B - 89 published spectra from 8 well-observed, 
nearby objects (see Table [5]). Additionally, the spectra 
were scaled to match the photometry. Given the absence 
of published uncertainties on these spectra, statistical 
errors were estimated from the pixel-to-pixel variation. 

The epochs used in this analysis are based on light 
curve estimates. In the case of Set A, the date of max- 
imum i?-band brightness was determined using prelimi- 
nary light curves. Therefore the epochs used to date the 
spectra were taken as the integer number of days since 
maximum light. The photometric data available for both 
set s were used to fit t he tem plate S-band light curve given 
by iGoldhaber et "aL (2001) and thus to obtain the values 
of the light curve decline parameters Am\^{B) and stretch 
(s) for each SN la. 

The two sets contain some peculiar SNe la, includ- 
ing the prototypes of t he two subclasses: SN 1991T and 
SN 1991bg. SN 1999 aa (|Garavini et al.ll2004h. SN 199 9aw 
(jStrolger et alJl2002h . and SN 1999bp (jFolatelli|2004j) are 
included in the 1991T-like subclass. All these SNe la 
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Table 4 SNe la from Set A (scp data: iFolateiil (|2004 )) 







Host 


Galaxy 


Rest-Frame 


SN 


z 


Galaxy 


Type a 


Spectral Epochs b 


1999aa 


0.0144 


NGC 2595 


SAB(rs)c 


-11, -7, -3, -1, 5, 6, 14, 19, 25, 28, 33, 40, 47, 51, 58 


1999ac 


0.0095 


NGC 6063 


Scd 


-15, -9, 0, 2, 8, 11, 16, 24, 28, 31, 33, 39, 42 


1999ao 


0.054 


Anon. 


S: 


5, 7, 10, 13, 18, 34, 40 


1999au 


0.124 


Anon. 


S: 


12, 17 


1999av 


0.05 


Anon. 


E/S0: 


2, 5, 9, 31 


1999aw 


0.038 


Anon. 


(?) 


3, 5, 9, 12, 16, 24, 31, 38 


1999be 


0.019 


Anon. 


(?) 


14, 19, 26, 33, 37, 44 


1999bk 


0.096 


Anon. 


E/S0: 


4, 7, 9 


1999bm 


0.143 


Anon. 


S: 


3, 6, 21 


1999bn 


0.129 


Anon. 


S: 


2, 14, 22 


1999bp 


0.077 


Anon. 


S: 


-2, 0, 1, 6, 17, 23 


1999bq 


0.149 


Anon. 


E/S0: 


3, 18 


1999by 


0.0021 


NGC 2841 


SA(r)b 


1, 6, 16, 27, 34 



a Hubble type of the host galaxy. An entry followed by a colon is a classification based on the host galaxy spectrum. 

rest are taken from NEE0 
b Rest-frame days since B-band maximum light. 



The 



Table 5 SNe la from Set B (Public data). 





2 


Host 


Galaxy 


Rest- Frame 




SN 




Galaxy 


Type a 


Spectral Epochs b 


References 


1981B 


0.0060 


NGC 4536 


SAB(rs)bc 


0, 17, 26, 29, 35, 49 


1 


1986G 


0.0018 


NGC 5128 


SO pec 


-7, -5, -1, ,3, 21, 28, 41, 44, 55 


2 


1989B 


0.0024 


NGC 3627 


SAB(s)b 


-7, -5, -3, -2, -1, 3, 5, 8, 9, 11, 16, 18, 19 


3 


1990N 


0.0034 


NGC 4639 


SAB(rs)bc 


-14, -7, 7, 14, 17, 38 


4, 5, 6 


1991bg 


0.0035 


NGC 4374 


El 


-2, 0, 13, 16, 23, 24, 30, 31, 44, 52 


7, 8 


1991T 


0.0058 


NGC 4527 


SAB(s)bc 


-11, -9, -8, -7, -6, -5, -3, 0, 10, 15, 23, 24, 42, 45 


9, 10, 11, 12 


1992A 


0.0063 


NGC 1380 


SAO 


-5, -1, 3, 7, 9, 11, 16, 17, 24, 28, 37, 46 


13 


1994D 


0.0015 


NGC 4526 


SAB(s) 


-10,-9,-8,-7,-5,-3, 2, 4, 6, 8, 11, 13, 14, 16, 18, 20, 25 


14, 15 



a Hubble type of the host galaxy from NED 5 . 

b Rest-frame days since -B- b and m aximu m light. 

i|l983h; (21 IPhillips et all i|l987l); (3) IWells et all 



Sourc es: (1) 
(5) iPhillips et al 



Branch et al 



(2) 

iMazzali et alj 



( 7) | FiIippenkp et al 



19921): (6) .. . .. 

(9) iFilippenko et all (ll992bT) ; (10) iRuiz-Lapuente et all dl992h: (1 1 ) IPhillips et al 
(13) iKirshner et ail l|l993h : (14) iMeikle et al] (|l996l ): (15) iPatat et al l (|l996h 




Leibundgut ct al 



Lcibundgut et al 
(12) iJefferv et al 



1991) 



1993) 



1992) 



present values of the decline rate parameter Amis(B) < 
1.0 (s > 1.1) and thus have a slow post-maximum de- 
cline in luminosity. At the other extre me, SN 1986G 
Phillips et alJll987t ICristiani et alJI 19921) and SN 1999by 



(jVinko et al.ll200ll : iGarnavich et alj|2004 belong to the 



1991bg-like subclass. These are fast-declining SNe la, with 
Aroi 5 (ff) > 1.70 (s < 0.80). The case of SN 1999ac 
( Garavini et al. 2005) is considered separately. This SN la 
has photometric and spectroscopic peculiarities that make 
it a unique object: its light curve shows a slow rise si milar 
to SN 1991T but a fast decline (jPhillips et alj|2002j) and 
its spectrum is similar to SN 1999aa. 



5 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 



4. Results 

We first perform a qualitative analysis (sections 14.11 and 
I4.2p and then turn to a more quantitative approach (sec- 
tions [O and ■ 



4.1. Absorption Velocities 

For normal SNe la, the magnitude of the blueshift of the 
Ca II H&K absorption minimum drops rapidly, from values 
around 22000 km/s before maximum light to 14000 km/s 
at maximum light. After maximum light, the decline in 
velocity flattens and decreases by about 4000 km/s in 50 
days. The mean trend for normal SNe la, from 10 days 
before maximum light to 40 days after maximum light, is 
shown in Fig. [4] together with the evolution in the veloc- 
ities of fast and slow-declining local SNe la. The shaded 
area represents the 1 
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sigma dispersion about the mean for normal SNe la. The 
trend and the dispersion hav e been computed from a large 
sample of nearby supernova (jGaravini et al.l l2004). 

The measured absorption velocity of the ejecta in 
under-luminous SNe la (i.e. SN 1999by and SN 1991bg 
plotted with the dotted and dashed lines, respectively) is 
approximately 1.5 sigma lower than the mean absorption 
velocity in normal SNe la. Therefore, the measurement of 
a low blueshift of the Ca n H&K absorption minimum 
cannot be used to identi fy a SNe la as under - luminous 
as already pointed out in iBlondin et al. (2006). Peculiar 
over-luminous SNe la such as SN 1991T and SN 1999aa 
show blueshifts of Ca n H&K absorption minimum that 
are within one sigma of those of normal SNe la. SN 1991T 
at one day before maximum light shows a blueshifts com- 
parable with that of SN 1999by, a fast-declining SN la. 

In Fig. [U the velocities of both the low-redshift and 
high-redshift SNe la are individually measured by per- 
forming an error-weighted non-linear least-squares fit to 
the entire line profile, manually selecting the end points 
limiting the wavelength region where to perform the fit. 
The line profile is modeled with a Gaussian plus a linear 
component. This method accurately reproduces the ab- 
sorption line profile a nd it has been s ucces sf ully used in 
previo us studies (e.g. Garavini et al.l ( 2005 ): Hook et al 



(|2005f l). All the SNe la in our data set were measured 
(Table El, with the exception of SN 2001gk, for which the 
line profile is incomplete, and SN 2001go at +29 days, for 
which the signal-to-noise is too low to correctly identify 
the absorption. The uncertainty in the redshifts is taken 
to be <j cz = 300 km/ s if the redshift was estimated from 
galaxy lines. For SN 2001ha and SN 2001hc, we could not 
identify host galaxy lines, so their redshifts were estimated 
by co mparing their spectr a with the spectra of nearby 
SNe la (jLidman et al.ll2005f ) . In these cases the uncertainty 
is increased to a z = 0.01 = a cz = 3000 km/s. This uncer- 
tainty dominates the statistical uncertainty from the fit. 
We note that the velocities of our high-redshift SNe la are 
consistent with the main trend of spectroscopically nor- 
mal local SNe la. In Fig [H open ci r cles i ndicate the mea- 
surements reported in Hoo k et alJ 1)20051 ) where a similar 
result was found in an independent sample of high redshift 
SNe la. 

4.2. Pseudo-Equivalent Widths 

Given the mean redshift of the high-redshift sample 
(z=0.49) and the typical S/N ratio of the spectra, we re- 
strict our analysis to the bluest and strongest features - 
#4 "Fe II 4800", #3 "Mg n 4300" and #1 "Ca II H&K". 

4.2.1. Estimating Systematic Errors 

Possible systematic errors arising from the choice of the fit- 
ting region at either side of the feature where the pseudo- 
continuum is fitted were accounted for by randomly shift- 
ing these regions (typically within a third of the region 




1-0x10* 



10 20 

Epoch referred to B mQ , (days) 



Fig. 4 The change in the blueshifts of the Ca n H&K 
absorption minimum with epoch for a sample of high- 
redshift SNe la (pre s ented in section [3j filled symbols, 
and from lHook et al. ( 2005 ) open circles) and a sample of 
low-redshift SNe la. The dashed and dotted lines indicate 
the values of extremely under-luminous SNe la SN 1999by 
(iGarnavich et alj|2004 ) and SN 1991bg (|Leibundgut et al 



19931 ) respectively. The solid line indicates the aver- 
age trend for Ca n H&K, which has been computed 
using a large data set of low-redshift normal SNe la 



(jGaravini et alJ 120041 ). The gray band shows the disper- 



sion (1 standard deviation) of the data about the average 
trend. 



size in each directional) and computing the weighted root- 
mean square deviation (rms) of the measured pseudo-EWs. 
This was the dominant source of uncertainties when the 
signal-to-noise ratio per resolution element was above ~10. 
Additionally, for lower signal-to-noise spectra (specifically 
in the case of our high redshift data-set) the central wave- 
length of the fitting region was randomly shifted accord- 
ing to a Gaussian distribution with er=10A. The change 
in the pseudo-EW is insignificant for high signal-to-noise 
ratio data. The standard deviations of these distributions 
were chosen in a conservative manner so as to take into 
account even large systematic effects. This source of error 
was added quadratically to the one given in Eqn. [2l 

It is known that pseudo-EWs can be affected by poor 
resolution a nd low signal-to-noise ratios (see, for example 
iGravl (|l992h ). Both effects were tested. Boxcar smoothing 
was used to decrease the resolution of the best-sampled 
spectra (~10 A/pixel) so that the range of resolutions in 
the present data set were tested. Due to the large intrinsic 
width of the broad SN la features, no significant change 
in the measured pseudo-EW was found. 

Signal-to-noise ratios ranged from about 5 to several 
hundred. When Gaussian noise was added to the best- 
quality spectra, in order to reproduce that quality range, 



6 Considering the uncertainty on the location of a region of 
size I to be Ijs/Vi ~ 1/3, as derived from the second moment 
of a uniform distribution. 
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Table 6 Magnitude of the blueshift of Ca n H&K absorption minimum. Measurement uncertainties are reported in 
parenthesis. 



SN 


day 


vel [km s i ] 


sn2001gy 


-7.5 (1) 


21520 (300) 


sn2002fd 


-7(2) 


14940 (300) 


sn2000fr 


-5.5 (1) 


17400 (300) 


sn2001gw 


-1 (2) 


15440 (300) 


sn2002gl 


0(2) 


13680 (300) 


sn2001ha 


0(2) 


16310 (3000) 


sn2001hc 


0(2) 


13110 (3000) 


sn2002gi 


2(2) 


13950 (300) 


sn2001gr 


2(2) 


15170 (300) 


sn2001gm 


5(2) 


17580 (300) 


sn2001go 


5.6 (1) 


14470 (300) 


sn2001go 


9.5 (1) 


13320 (300) 



no significant bias was detected in the resulting pseudo- 

EWS. 



Additionally, the effect of reddening was tested by ar- 
tificially adding up to Eb-v = 0.32 magnitudes of red- 
dening (corresponding to Ay = 1 mag, with Ry = 3.1) 
following the law given in lCardelli et al. (1989). This pro- 
duced no significant change in the resulting pseudo-EWs. 
This is expected since the pseudo-EWs are defined over a 
limited wavelength range. 



Table 7 The fractional decrease in the pseudo-EW corre- 
sponding to a 10% increase in the amount of contamina- 
tion from the host. 



Further systematic effects could arise from host-galaxy 
light that has not been perfectly removed. The effect of 
additional signal underlying the SN la spectrum would be 
to lower the pseudo-EW. The spectra from Set B in the 
present sample correspond to very bright, nearby SNe la, 
for which SN la and host-galaxy spectra can be resolved. 
For more distant SNe la, from Set A and the high red- 
shift data set, the light from the host can contribute up 
to 50% of the light in the extracted spectra. We have 
tested how errors in estimating the amount of host galaxy 
light can affect the pseudo-EWs of a typical near-maximum 
light SNe la spectrum. Template galaxy spectra of Hubble 
types E and Sc were added to the SN la spectrum in or- 
der to simulate contamination levels ranging up to 50% of 
the total integrated flux between 3500 and 9000 A. The 
pseudo-EWs of the features were then measured on every 
spectrum. The relative decrease in the pseudo-EW with 
increasing contamination levels was found to be approxi- 
mately linear. Table [7] summarizes these results by giving 
the relative decrease in the pseudo-EW per each 10% of un- 
accounted contamination on the total integrated flux for 
both galaxy types. Since SNe la near maximum light are 
generally bluer than their hosts, errors in estimating the 
amount of host galaxy contamination leads to larger errors 
in the pseudo-EWs for features at redder wavelengths. For 
early type galaxies, because of the presence of the Balmer 
break around 4000A, the effect on the pseudo-EW of "Ca II 
H&K" is less than 10% even for 50% contamination. 



Feature 


Mnemonic 


Host Type 


ID 


Label 




Sc a 


1 


"Ca n H&K" 


0.019 


0.080 


2 


"Si ii 4000" 


0.048 


0.074 


3 


"Mg ii 4300" 


0.037 


0.073 


4 


"Fe ii 4800" 


0.070 


0.074 


5 


"S ii W" 


0.112 


0.074 


6 


"Si ii 5800" 


0.128 


0.066 


7 


"Si ii 6150" 


0.103 


0.084 


8 


"Ca ii IR" 


0.155 


0.114 



Fractional decrease in pseudo-EW. 



4.2.2. "Fe ii 4800" (#4). 

The evolution in the pseudo-EW of feature #4, "Fe n 4800", 
in nearby supernovae is shown in Fig. [5l The number of 
data points enable us to compute a mean trend, which is 
shown as the solid line in Fig. [5l for normal SNe la only. 
The curve was built in the range — lOd < epoch < 50d 
by dividing the epochs into 5 bins, calculating a weighted 
average of the pseudo-EW in each bin, and finally tracing 
the spline function through those points as a general in- 
dication of the followed time evolution. The coordinates 
(epoch, pseudo-EW), in units of days and A, of the 5 points 
defining the curve are: (-5,134); (5,181); (15,267); (25,339); 
(35,356). 

The pseudo-EW of "Fe II 4800" monotonically increases 
with phase from before maximum light to about 30 days 
after maximum light. This is due to the increasing optical 
depth of Fe II lines from around maximum light onward 
and to the subsequent overlapping of several Fe II lines 
from around 15 days after maximum light (see Fig. []}. 

The evolution of pseudo-EWs is similar for all la sub- 
types, but offset from the mean curve. In general, 1991bg- 
like SNe la lie above the curve and 1991T-like SNe la lie 
below it but there is no firm correlation with Amis. This 
is summarized in Table |8] which shows the distribution 
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of the three SNe la subtypes with respect to the average 
curve. 

We note that for normal type la SNe the average trend 
in the first data bin is based on four data points from 
SN 1994D and one each from SN 1989B and SN 1990N. 
The mean is then biased toward the average pseudo-EW 
value of SN 1994D between -10 and -6 days. This is higher 
than the measurement obtained on the other SNe la, thus, 
the pre- maximum average trend is biased toward the trend 
observed in SN 1994D. More pseudo-EW measurements are 
needed to firmly establish the pre-maximum evolution in 
normal SNe la. 

In Fig [5] (lower panel) the "Fe II 4800" pseudo-EWs 
of the high redshift sample are shown. All SNe la (with 
the exception of SN 2001gu, SN 2001gw, SN 2001gy and 
SN 2002gi, for which the absorption feature was not easily 
identifiable) were found to lie within the 95% probability 
distribution of low-redshift supernovae indicated by the 
gray filled area. 

4.2.3. "Mg ii 4300" (#3) 

Various ions correspond to feature #3. These include 
Mg II, Co II, Fe II, Fe III, and Si in for spectroscopically 
normal and 1991T-like SNe la. In the case of 1991bg-like 
SNe la, the region is dominated by strong lines of Ti II 
(|Filippenko et al.lll992al : iMazzali et al.lll997t ). The evolu- 
tion of the pseudo-equivalent width of this feature is dif- 
ferent than that of "Fe II 4800" as can be seen in Fig. El 

The pseudo-EW of the feature increases dramatically 
over a short period of time as it mergers with the neigh- 
boring "Si II 4000" feature (feature #3 in figure [5]). Before 
and after this increase, the pseudo-EW of this feature is 
approximately constant. The phase at which this increase 
takes place, t{, r , is highly dependent on the SN la sub- 
type. For 1991bg-like SNe la it seems to occur as early 
as 5 days before maximum light (the earliest spectrum of 
a 1991bg-like SN la in our sample), while normal SNe la 
show this behavior around one week after maximum light, 
and 1991T-like objects show it later than day +10. Thus, 
the evolution of the pseudo-EW of the "Mg n 4300" feature 
can be used to discriminate between different Type la sub- 
types. 

We describe the average evolution of this feature with 
the function: 



f(0,t) 



A 



1 



(3) 



with parameters 6 — (A, B, tt, r , t). Parameters (A, B) 
are simple constants, t is the SN phase, and r is an e- 
folding time scale. The solid line in Fig. [H] shows the aver- 
age curve for normal SNe la between days —10 and +30. 
Other SN la subtypes show a parallel trend. Table O lists 
the dispersion of normal SNe la around the average curve 
and the deviations of 1991T-like and 1991bg-like SNe la 
from the same curve. The parameterization introduced in 
Eq. [3] allows one to quantify the sudden change in the 



□ 1 99 1 bg — like 
• normals 
o 1 99 1 T— like 
A SN 1999ac 
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Fig. 5 "Fe II 4800" 

Upper Panel: Measured pseudo-equivalent width corre- 
sponding to the "Fe II 4800" feature (#4). SN 1991bg-like 
objects are marked with open squares, 1991T-like SNe la 
with open circles, normal SNe la with filled circles and 
SN 1999ac with triangles. Error-bars include the error de- 
scribed by Eqn. [2] as well as systematic uncertainties aris- 
ing from the pseudo-continuum fit. The solid line shows 
a cubic spline function used to represent the average evo- 
lution of normal SNe la between days —10 and +50. In 
general, 1991bg-like SNe la lie above the average curve 
whereas 1991T-like SNe la lie below it. Lower Panel: A 
comparison between the 'Fe II 4800' pseudo-EWs in low- 
and high-redshift SNe la. High-redshift supernova are in- 
dicated by large filled symbols. The gray filled area repre- 
sents the 95% probability region for normal low-redshift 
SNe la. See text for details. 

pseudo-EW in terms of the parameter t^ r , which can then 
be used to classify type-la supernovae into the three main 
sub-types. In Fig. [5] (upper panel) we note that, for nor- 
mal type la SNe before maximum light, the pseudo-EW 
measurements tend to cluster around high and low val- 
ues. The high values of pseudo-EW are those measured 
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Table 8 Dispersion of "Fe n 4800" pseudo-EW for the three SNe la subtypes. See text for details. 



Epoch 




Normal 






1991T-like 


1991bg-like 


bin 


n 


(pseudo-Ew) 
(A) 


6 

(A) 


n 


^pscudo-EW 

(A) 


n 


•^pscudo-EW 

(A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


[-10,-6] 


6 


132 ±11 


28 


5 


-29 ±3 







[-5,-1] 


6 


114 ± 12 


37 


5 


-36 ± 5 


2 


16 ±4 


[0,5] 


12 


164 ± 10 


28 


5 


-62 ±9 


4 


24 ± 7 


[6, 11] 


15 


200 ±8 


25 


4 


-58 ±15 


1 


113 a 


[12, 18] 


15 


265 ±9 


25 


5 


-47 ± 12 


1 


139 a 


[19,27] 


10 


316 ± 12 


30 


5 


-34 ± 18 


3 


58 ±8 


[28, 40] 


9 


353 ± 13 


36 


5 


3±9 


3 


37 ±5 



Columns: (1) Range of epoch bins in days; (2) Number of points from normal SNe la; (3) Average pseudo-EW for normal 
SNe la; (4) Dispersion (rms) of normal SNe la around the average curv^ll; (5) Number of points from 1991T-like SNe la; 
(6) Mean deviation of 1991T-like SNe la from the average curve. Uncertainties do not include the computed dispersion of 
normal SNe la around the curve (8); (7) Number of points from 1991bg-like SNe la; (8) Mean deviation of 1991bg-like SNe la 
from the average curve. Uncertainties do not include the computed dispersion of normal SNe la around the curve (5). 
a Only one measurement, thus no uncertainty is given. 



on SN 1989B (two data points in the first epoch bin and 
three in the second one) and on SN 1990N (one data point 
in each bin). Instead, SN 1994D shows low pseudo-EWs 
(two data points in each bin before maximum light), with 
SN 1992A pseudo-EW having intermediate values (in both 
pre-maximum bins). The average trend is then greatly af- 
fected by SN 1994D which, as for feature #4, seems to be 
an outlier. Because of the low statistics the pre-maximum 
trend is only indicative. An analysis including more nor- 
mal type la SN will be needed to assess the average trend 
and to investigate whether the apparent clustering in the 
distribution is a distinctive characteristic of type la SNe 
disclosing interesting physics. 

From Fig. [5] (upper panel), it is evident that the evolu- 
tionary behavior of the "Mg n 4300" pseudo-EW is different 
for normal, 91T-like, and 91bg-like SNe la, and therefore 
could correlate also with photometric properties of the SN. 
To quantify this correlation, the functional model given by 
Eq.[3]was used to fit the parameter t\, r for each of SN la in- 
dividually. This parameter is related to the phase at which 
this feature suddenly becomes stronger. Table [TU1 lists the 
values of tf, r for 11 SNe la in our sample. In the cases of SN 
1991bg and SN 1999by, the epochs of their earliest spec- 
tra are used as upper limits for tb r , assuming they followed 
the same evolutionary pattern found for the other objects. 
The values of Atoi5(B) were re-fitted for this analysis 
to gain consistency in the results. Fig. [7] shows the cor- 
relation between Ami^(B) and ti, r . The break occurs at 
later epochs for SNe la with slower declining lightcurves. 
A least-squares fit yields: 



Ami 5 (B) = 1.869(±0.052) - 0.070(±0.005)i 6r , (4) 



7 In general we find S ~ s/N ■ <5p S eudo-EW> ^ut n °t e that 
there are slight deviations due to the use of different epoch 
bins for this table and in building the average trend. 



where 1.7 < i& r < 16.9 is given in days since maximum 
light. The fit was done excluding the two upper limits. The 
resulting dispersion in Amis(i?) around the fitted line is 
0.08 mag, which is comparable to the measurement un- 
certainties. Since this spectral feature is covered by the 
standard B filter, the correlation is not surprising. The 
sudden strengthening of the "Mg II 4300" absorption ap- 
pears to affect the shape of the lightcurve. 

An additional point to note is that this feature can 
be observed with optical spectrographs up to redshifts 
of ~ 1, and the technique can be used up to redshifts 
> 1.7 with future SN-related cosmology experiments (e.g 
ADEPT, DESTINY or SNAP) which will include near- 
infrared spectrographs. However, since several spectro- 
scopic observations are needed to obtain the value of t{, r , 
the use of this parameter to estimate Mg, lax would be lim- 
ited to SNe la with rather extensive spectroscopic follow- 
up. 

FigE] (lower panel) shows that the "Mg II 4300" pseudo- 
ews of the high-redshift supernovae (all except SN 2002gi 
are measured) are consistent with the trend defined by 
the 95% probability distribution of low-redshift normal 
supernova indicated by the gray filled area. Moreover, 
SN 2001go, for which measurements at three epochs are 
available, shows the sudden increase at around one week 
after maximum, as do most low-redshift normal SNe la. 

4.2.4. "Ca II H&K" (#1) 

The prominent absorption trough at 3800A is attributed 
to the Ca n H&K lines with contributions from Si n 
A3858 and lines from iron-peak elements. Fig. [5] shows 
the change in the pseudo-EW of this feature with phase. 
As expected for SN 1991T-like objects, this absorption is 
particularly weak (low pseudo-Ew), especially in the pre- 
maximum spectra. The intrinsic dispersion in the pseudo- 
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Table 9 Dispersion of "Mg n 4300" pseudo-EW for the three SN la subtypes. See text for details. 



Epoch 




Normal 






1991T-like 


1991bg-like 


bin 


n 


(pseudo-Ew) 


6 


n 


^pscudo-EW 


n 


^pscudo-EW 






(A) 


(A) 




(A) 




(A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


[-10,-5] 


6 


90 ±7 


16 


6 


15 ±4 







[-4,0] 


6 


98 ±5 


14 


5 


-5 ±3 


4 


139 ± 22 


[1,5] 


9 


118 ± 12 


28 


4 


-25 ±12 


1 


206 a 


[6, 10] 


10 


184 ± 14 


26 


4 


-77 ± 12 







[11, 17] 


15 


251 ± 10 


33 


5 


-69 ±23 


2 


48 ± 2 


[18, 30] 


15 


257 ±7 


27 


6 


8±6 


4 


21 ±20 



See column explanations in Table [8] 
a Only one measurement, thus no uncertainty is given. 



Table 10 AmiQ(B) and tb r measured on the nearby supernova sample. 



SN 


Am 15 (B) a 


i b 

b\y T 


1986G 


1.73(05) 


1.7 


1989B 


1.32(07) 


7.1 


1990N 


1.03(10) 


9.2 


1991T 


0.96(07) 


12.1 


1991bg 


1.93(05) 


< C 


1992A 


1.45(04) 


7.0 


1994D 


1.46(10) 


7.0 


1999aa 


0.95(08) 


13.6 


1999aw 


0.75(10) 


16.4 


1999bp 


0.70(06) 


16.9 


1999by 


1.90(05) 


< l c 



a Uncertainties given between brackets, in units of 0.01 mag. 

b Spectral parameter related to "Mg II 4300" (Eqn. [3}. Assumed uncertainty in tb r : 0.5d. Measured in days since maximum 
light. 

c Upper limit corresponding to the earliest data point. 



EWs is greater than in the case of features #3 (Mg II 4300) 
and #4 (Fe n 4800). 

This is primarily due to the strength of the Ca n H&K 
line when compared with Mg II and Fe n . The evolution 
of the pseudo-EW with phase is specific to each object 
both qualitatively and quantitatively. The dispersion on 
pseudo-EW is greater before maximum light, where some 
objects show an increase and others a decrease. After max- 
imum light, there is normally a slow decrease, except for 
the 1991bg-like objects in which the pseudo-EW is rela- 
tively constant. The spread of pseudo-EWs of this feature 
is the largest of all the features analyzed in this work. 
This is due, in part, to the possible presence of high veloc- 
ity components and to the overlap of several lines that lie 
in the wide wave length interval spann ed by this feature. 
As pointed out in Ma zzali et al. ( 20051 ) the position in ve- 
locity space of such components differ in different SNe la 
making the feature broader or narrower and thus affecting 
the measured pseudo-EW. 

The two points with pseudo-EW > 170 A after day 
belong to SN 1999bm. Their large pseudo-EW values are 
due to an unusually br o ad "C a n H&K" feature. As pointed 
out in Gerardv et al. ( 20041 ) this phenomenon might be 



caused either by the presence of a strong, high-velocity 
Ca II component or by a strong Si II 3858 line. Disentangle 
between these two hypotheses is difficult. The presence of 
a high-velocity Ca II component could be better deter- 
mined in the Ca n IR triplet region, because of the lack of 
contamination from other ions. However, in our spectra of 
SN 1999bm the signal-to- noise around 8000 A is two low 
to identify such component. It is impossible to exclude 
an high-velocity component based only on Ca II H&K. 
The strength of the Ca II H&K line could result in a sin- 
gle minimum absorption line even when an high-velocity 
components is present and detected in the Ca II IR triplet. 

The pseudo-EW of "Ca II H&K" feature of low- and 
high-redshift supernova are shown in Fig. [5] For this fea- 
ture, as described above, it is not possible to identify 
a distinctive trend for SN la subgroups as a function 
of time. Before twenty days past maximum, normal and 
under-luminous low-redshift SNe la populate the region 
that spans from pseudo-E\v60 A to pseudo-E\vl40 A. In 
the lower panel of Fig. [SJ the gray filled area represents 
the 95% probability region for normal (non-outliers) and 
under-luminous low-redshift SNe la. Some outliers (deter- 
mined by a 3-sigma clipping algorithm) are found among 
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Fig. 6 "Mg II 4300" 

Upper Panel: Measured pseudo-equivalent widths corre- 
sponding to "Mg ii 4300" (#3). SN 1991bg-like objects are 
marked with open squares, 1991T-like SNe la with open 
circles, normal SNe la with filled circles and SN 1999ac 
with triangles. Error-bars include the error described by 
Eqn.[5]as well as systematic uncertainties arising from the 
pseudo-continuum fit. The solid line represents the aver- 
age behavior of normal SNe la, as described in Eqn.[3]with 
0=(172.2, 89.4, 7.5, 2.5). In general, 1991bg-like SNe la 
lie above the average curve whereas 1991T-like SNe la 
lie below it. Lower Panel: A comparison between the 
'Mg II 4300' pseudo-EWs in low- and high-redshift SNe la. 
High-rcdshift supernova are indicated by large filled sym- 
bols. The gray filled area represents the 95% probabil- 
ity region for normal low-redshift SNe la. Peculiar under- 
luminous nearby SNe la are shown separately for compar- 
ison. See text for details. 

normal SNe la and are indicated with small diamond 
symbols. Before maximum light, peculiar 91T-like objects 
show systematically low values, as indicated by the small 
square symbols. The high-redshift supernovae (all except 
SN 2002gi are plotted) do not show significant deviations 




5 10 

t br (days since B maximum) 



Fig. 7 The "Mg n 4300" pseudo-EW break parameter U r 
(Eqn.[3]). Amis(B) versus tt, r . The two points marked with 
arrows correspond to upper limits of tb r . The straight line 
is least-squares fit to the data, excluding the upper limits. 



with respect to the low-redshift sample shown in the plot, 
and SN 2002fd falls on the 91T-like trend as expected. 

4.3. The "Ca II H&K" pseudo-EW of SN 2002fd 

Prior to maximum light, the pseudo-EW of Ca II H&K 
can be used to separate 91T-like SN la from normal 
SNe la, (see Fig. If the identification of SN 2002fd 
as a peculiar object is correct, we expect the pseudo-EW 
to be lower than that of normal SNe la. The average 
pseudo-EW prior to maximum light in normal SNe la is 
<pseudo-EW> = 114.1 and the scatter around the mean 
value is <7<p SC udo-Ew>= 14.2 For peculiar 91T-like SNe la 
we find <pseudo-EW> = 68.7 and cr< pS cudo-Ew> =6.1 The 
value measured for SN 2002fd, (pseudo-EW = 73.6 ± 2.9), 
is consistent — within one standard deviation — with that 
found for 91T-like SN la, and inconsistent (at more than 
3 standard deviations) with normal SNe la. 

Note, that the plotted errors bars in the pseudo-EWs 
of the high-redshift SNe la include both statistical uncer- 
tainties, from the measurement, and systematic uncertain- 
ties from residual host galaxy contamination but do not 
include systematic uncertainties from possible misidentifi- 
cation of the maxima around the absorption feature (i.e. 
misidentification of the fitting regions). Table [Til reports 
the measured pseudo-EWs. 

4.4. Statistical comparison 

In this section, a statistical comparison of the low- and 
high-redshift SNe la is performed using the spectral in- 
dicators described in section [5] The mean trends in the 
pseudo-EWs of the 'Fe n 4800' and 'Mg n 4300' features 
identified for normal - low-redshift - SNe la can be used to 
test whether or not high-redshift supernovae pseudo-EWs 
follow the same trends. 
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Table 11 Measurements of the pseudo-EWs of Ca n H&K, Mg II and Fe II. Measurements uncertainties are reported 
in parenthesis. a 



SN 


epoch 


Ca n H&K pseudo-EW 


Mg II pseudo-EW 


Fe II pseudo-EW 




[days] 


[A] 


[A] 




[A] 


sn2001gy 


-7.5 (1) 


176 (8/21) 


111 (12/31) 






sn2002fd 


-7(2) 


73.6 (3/14) 


101 (8/11) 


86 


(15/23) 


sn2000fr 


-5.5 (1) 


122 (13/14) 


103 (6/13) 


91 


(13/23) 


sn2001gw 


-1 (2) 


130 (9/13) 


125 (17/22) 






sn2002gl 


0(2) 


124 (14/15) 


121 (12/18) 


168 


(20/25) 


sn2001ha 


0(2) 


188 (19/22) 


132 (9/24) 


216 


(33/62) 


sn2001hc 


0(2) 


112 (15/15) 


116 (11/18) 


140 


(15/23) 


sn2001gr 


2(2) 


94.6 (10/22) 


140 (16/34) 


141 


(31/69) 


sn2001gm 


5(2) 


128 (23/26) 


123 (17/40) 


169 


(34/84) 


sn2001go 


5.6 (1) 


120 (11/28) 


93 (11/24) 


191 


(28/41) 


sn2002gk 


6(2) 


91 (6/11) 


161 (14/16) 


157 


(16/17) 


sn2001go 


9.5 (1) 


112 (6/13) 


181 (14/24) 


174 


(20/27) 


sn2001go 


29.5 (1) 


228 (22/48) 


315 (41/52) 







a The first value includes statistical and systematic uncertainties related to the level of host galaxy contamination, the second 
value (added in quadrature) also includes possible systematic uncertainties due to misidentification of the fitting regions, (see 
sec. [2] for details). 



The results of a set of % 2 tests are shown in Table 
[T2l The intrinsic dispersion around the fitted mean trends 
for normal low-z supernovae was added in quadrature to 
the statistical and systematic uncertainties to perform the 
test. The uncertainty in the SN la phase was propagated 
according to the pseudo-EW model for normal low redshift 
SNe la and added in quadrature to the measurements er- 
ror on the pseudo-EWs. The possible systematic uncertain- 
ties due to misidentification of the maxima was included 
was also taken into consideration. We note that the fit- 
ting region uncertainties — included in the results shown 
in Table [T2l — should be considered as upper limits to the 
possible systematic uncertainty due to the pseudo-EW's 
measurement technique (see section [5] for details). 

The hypothesis that the pseudo-EW measured on our 
high-redshift supernovae follow the same behavior with 
lightcurve phase as those measured on low redshift normal 
supernovae is statistically confirmed. Moreover, the hy- 
pothesis that pseudo-EWs measured on our high-redshift 
supernovae are consistent with those of under- luminous 
SN 1991bg-like low redshift SNe la is rejected. 

5. Summary and Conclusions 

Spectroscopic data of 12 high-redshift supernovae, in the 
redshift interval z= 0.212 to 0.912, were analyzed and a 
qualitative classification scheme was presented. Based on 
this classification scheme, all of our high-redshift SNe la 
were classified as normal SNe la, except for SN 2002fd 
(z=0.27), which is similar to SN 1999aa, a peculiar 91T- 
like SN la. We also find, based on spectral properties 
alone, that none of the supernovae studied in this work 
are under- luminous. This is not unexpected because of 
the bias against selecting suc h objects in mag nitude lim- 
ited surveys (see for example iLi et al.l (|2001al )). 



A quantitative comparison between low and high- 
redshift SNe la by means of spectral indicators has been 
presented. The velocities of the minimum of Ca n H&K 
feature of high-redshift SNe la, were compared to those of 
low-redshift SNe la with the aim of uncovering differences. 
No systematic differences could be found. 

Using a low-redshift SN la sample, we study the evo- 
lution in the pseudo-EW of the strongest spectral features 
as a function of phase, and find that the pseudo-EW of dif- 
ferent SNe la sub-types (normal, 91T-like and 91bg-like) 
evolve differently for two of the three features studied (i.e. 
"Fe II 4800" and "Mg n 4300"). In the case of "Fe n 4800", 
the different SNe la sub-types follow similar trends, but 
offset from the average evolution. 91bg-like SNe la show 
higher pseudo-EWs and 91T-likc SNe la lower pseudo-EWs 
than normal SNe la. The pseudo-EW of "Mg n 4300" is 
characterized by a sudden break around maximum light. 
The epoch at which the break occurs correlates with the 
photometric properties of the SNe la. We find that in 
91bg-like SNe la the break occurs earlier in phase with 
respect to normal SNe la while in 9 IT- like SNe la the 
same occurs later. 

The pseudo-equivalent widths of "Fe n 4800", "Mg n 
4300" and "Ca n H&K" in high-redshift SNe la are found 
to follow the same trends with epoch as those observed 
in normal low-redshift SNe la. Furthermore, the pseudo- 
equivalent widths of "Fe II 4800" and "Mg n 4300" in high- 
redshift SNe la are found to be statistically consistent with 
the pseudo-equivalent widths observed in low-redshift nor- 
mal SNe la. 

The pseudo-EWs of Ca n H&K in the spectrum 
of SN 2002fd are consistent with those observed in 
SN 1991T/SN1999aa-like objects in the local Universe, 
quantitatively confirming the sub-type identification of 
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Table 12 A statistical comparison of the 'Fe n 4800' and 'Mg n 4300' pseudo-equivalent widths of high and low-redshift 
SNe la. 



Feature 


Xnorm 


Pnorm 


'2 

X91bg-likc 


P91bg-likc 


Xnorm 


Pnorm 


•2 

X91bg-like 


P91bg-liko 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Fe n 


13.0(11) 


0.29 


29.6(11) 


0.002 


6.7(11) 


0.8 


16.8(11) 


0.1 


Mg ii 


11.2(13) 


0.59 


228.4(13) 


0.00 


7.3(13) 


0.9 


162.4(13) 


0.00 



Columns 2, 3, 6 and 7 report the x 2 correspond to the comparison of the pseudo-EWs 'Fe II 4800' and 'Mg II 4300' of high 
redshift SNe la with the fits for normal supernovae. Columns 4, 5, 8 and 9 report the comparison with under-luminous SNe la 
(e.g. SN 1986G and SN 1991bg). Columns 2, 4, 6, and 8 reports the \ 2 obtained with the number of degrees of freedom 
in parenthesis. Columns 3, 5, 7, and 9 report the probability to obtain a x 2 value greater than that obtained. Columns 2 
to 6 report the results including both statistical uncertainties and systematic uncertainties due to host galaxy subtraction. 
Columns 7 to 9 report the results when the systematic uncertainty due to possible misidentification of fitting regions is also 
included. 



this SN la. This feature can be used to identify 91T-like 
objects at high-redshift. 

The number of higher-redshift supernova used in this 
case study is small when compared to the numbers of 
SNe la that arc now being observed at similar redshifts. 
Both the ESSENCE and SNLS projects will result in sam- 
ples of SNe la with similar spectral quality; however, the 
samples will be many times larger. The case study pre- 
sented here offers a simple method to analyze the spectra 
observed in these surveys and look for systematic differ- 
ences between SNe la at different redshift. 
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Fig. 8 "Ca ii H&K" 

Upper Panel: Measured pseudo-equivalent widths corre- 
sponding to "Ca ii H&K" (#1). SN 1991bg-like objects are 
marked with open squares, 1991T-like SNe la with open 
circles, normal SNe la with filled circles and SN 1999ac 
with triangles. Error-bars include the error described by 
Eqn.[5]as well as systematic uncertainties arising from the 
pseudo-continuum fit. The vertical dashed line indicates 
the upper epoch limit in the lower panel. Lower Panel: A 
comparison between the 'Ca II H&K' pseudo-EWs in low- 
and high-redshift SNe la. High-redshift supernova are in- 
dicated by large filled symbols. The gray filled area repre- 
sents the 95% probability region for normal (non-outliers) 
and the under-luminous low-redshift SNe la SN 1999by,SN 
1991bg and SN 1986G. Outliers normal SNe la are indi- 
cated with the diamond symbols. Slow declining SNe la 
are indicated by the open square symbols. 
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